Quantum Shock Waves - 
the case for non-linear effects in dynamics of electronic liquids 
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Using the Calogero model as an example, we show that the transport in interacting non-dissipative 
electronic systems is essentially non-linear. Non-linear effects are due to the curvature of the elec- 
tronic spectrum near the Fermi energy. As is typical for non-linear systems, propagating wave 
packets are unstable. At finite time shock wave singularities develop, the wave packet collapses, and 
oscillatory features arise. They evolve into regularly structured localized pulses carrying a fraction- 
ally quantized charge - soliton trains. We briefly discuss perspectives of observation of Quantum 
Shock Waves in edge states of Fractional Quantum Hall Effect and a direct measurement of the 
fractional charge. 
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1. It is commonly assumed that a small and smooth 
perturbation of electronic systems remains small and 
smooth in the course of evolution. This assumption jus- 
tifies the linear response theory to electronic transport. 
Although it is true for dissipative electronic systems, like 
a metal, it fails for dissipationless quantum fluids. 

In this paper we argue that in such electronic systems, 
commonly known in one spatial dimension, the trans- 
port is essentially non-linear and is subject to an unsta- 
ble singular behavior. An arbitrarily small and smooth 
disturbance of electronic density and momentum (a wave 
packet) will inevitably develop to shot-noise-like oscilla- 
tory features, progressively evolving to regular localized 
pulses, carrying quantized charge - a soliton train. This 
process is a dissipationless counterpart of a shock wave, 
where an initially smooth profile of density tends to over- 
hang and collapse at some finite time. A typical behavior 
is shown in Fig.^ We will demonstrate this phenomenon 
on the example of the Calogero model, but believe that a 
shock wave instability and the main features of the sub- 
sequent dynamics are universal for dissipationless inter- 
acting electronic systems. They are insensitive to details 
of interactions and to the initial profile of the wave. 

We note that, for some physical systems, the time of 
the onset of the shock wave may be longer than the time 
needed for other physical mechanisms to dissipate the 
disturbance. However, our estimates indicate that the 
observation of shock waves at confined states of a quan- 
tum well heterostructures (the edge states of Fractional 
Quantum Hall Effect (FQHE) ) is possible. In addition, 
a dynamic experiment may allow direct observation of 
the fractional charge of "fly-out" excitations. 

2. Calogero-Sutherland model defined on a circle by 
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occupies a special place in ID quantum physics. The sin- 
gular interaction is so strong that it changes the statistics 
of particles from 1 (Fermi-statistics) to A, while excita- 
tions carry statistics of 1/A (for a review see e.g., QQ)- 
The chiral sector of Calogero model, therefore, describes 
the edge of the FQHE with a fractional filling 1/A. Here 
A is a dimensionless coupling constant g = — A(A — 1). 
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FIG. 1: Numerical solution of BO equation JHj- A: The initial 
Lorentzian profile containing 7 solitons is shown at time t = 
0.9i c . B: The wave after a shock is shown together with the 
overhanging solution of Riemann equation JSJ (red, dashed) 
at t — 2.7t c . Vertical lines mark the trailing and the leading 
edges. The moduli v,v,K of the 1-phase solution @ are 
assigned to different branches of the multi- valued solution of 
©. D: The solutions of BO and Riemann's equations are 
shown at t = llt c . C: Whitham modulated wave represents 
an approximate soliton train behavior for an initial Lorentzian 
shape of a larger area. 
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3. Evolution of a semiclassical state. Consider an 
evolution of an initial state - a wave-packet |\&), of a 
"large" size I » kp 1 as on Fig. The state is as- 
sumed to be semiclassical, i.e., it involves many parti- 
cles, and is such that its Wigner's function W(x,p) — 
^\ e i( Px +P x )\-qi) has a well-defined classical limit. This 
means that W(x,p) exponentially vanishes outside of a 
certain star-like support in phase space. Semiclassical 
hydrodynamic modes can be viewed as area-preserving 
waves of the boundary of the support. We study the 
time evolution of the density p{x, t) = JW(x,p, t)dp and 
current j{x,t) = J pW(x,p,t)dp, or velocity v{x,t) de- 
fined as j = mpv of the wave-packet. The state can be 
created by some classical instrument acting on a system. 

First, we argue that in interacting ID system like 
the evolution of a semiclassical state is determined only 
by hydrodynamic variables - the density p(x,t) and ve- 
locity v(x,t). In other words, the semiclassical density 
and velocity obey a set of closed evolution equations 
of hydrodynamic type. Things arc different for non- 
interacting particles A = 1. There, a "gradient catastro- 
phe" also occurs at a finite time, but further evolution 
has an essentially quantum nature. 

4- Quantum hydrodynamics. To study the dynam- 
ics of and similar problems, it is productive to rep- 
resent a quantum many body problem as quantum hy- 
drodynamics, i.e., rewrite it solely in terms of the den- 
sity p(x) — J2iH x ~ and the velocity operators 
[p(x), v(y)] — —i^S'(x-y) 4]. In this approach [5| excita- 
tions appear as solitons of nonlinear fields p{x) and v(x). 
Calogero model has been successfully written in terms 
of hydrodynamic variables, thus extending the familiar 
"bosonization" of fermionic models in ID 13, EJ 0, 0] • 
Below we remind the major facts following [3j . 

Equations of motion of JQ) are equivalent to equations 
of quantum hydrodynamics 

continuity equation : p + 'V(pv) =0, (2) 

Euler equation : v + V(^ + w) = 0, (3) 

where w = (K/2ir) 2 5(pe)/5p is the quantum enthalpy 
and e is the energy per particle in units of \k/2it 

la 2 1 
e= -W 2 + - (Vlnp) 2 + -™Vp ff , (4) 

and pn = ^p{x') cot x) is a Hilbert transform 

of the density. The first term in (@J is present already for 
free Fermi gas (A = 1) - the Fermi energy. The last two 
terms absent at A = 1 reflect the interaction. 

5. Sound waves, nonlinear effects and dispersion. In 
the limit of small deviations of density Sp = p — po <C po, 
and velocity Sv = v — vq <C vq from their mean values 
po and vq, eqs. (12131) may be linearized. If in addition 
gradients are small a-^p <C po so that the two last terms 
in can be dropped, one obtains the familiar "linear 



bosonization" , where density and velocity waves move 
to the right (left) with Fermi velocities vf = «o ± §Po 
without dispersion. Introducing left and right moving 
fields u = Sv ± ^Sp the linear hydrodynamics reads u + 
j) s Vm = 0, where v s — vf are velocities of sound. 

Shifting to a Galilean frame moving with velocity v s , 
one recovers non-linear effects: 

it + uVu = 0. (5) 

This is is the quantum Riemann equation for ID com- 
pressible hydrodynamics [lj}. Its follows from the 
Galilean invariance and may be traced to a curvature 
of the particle spectrum. The Riemann equation states 
that the wave u{x 1 1) moves with a velocity which is itself 
given by u(x,t). Contrary to the case of a linear spec- 
trum, when a packet moves with the velocity of sound, 
and does not change its shape, an arbitrarily small disper- 
sion deforms the packet pushing its denser part forward. 

Eq. © is equivalent to ill'l.'il) under assumption of small 
gradients a\/Sp -C Sp 2 when the two last terms in @ are 
small. However this assumption almost never holds. It 
fails for any small and smooth initial condition with a 
"bright side" when u(x)Vu(x) < (Fig. E}\). This is 
seen from the implicit solution of JSJ due to Riemann: 
u = f(x — ut), where f(x) = u(x,0) is the initial state. 
The bright side of the wave packet steepens and finally 
overturns at some time, t r , where the formal solution 
becomes multi- valued jlClj . and un-physical Fig. 03-D. 

Hydrodynamics of the electronic system can be un- 
derstood as a hydrodynamics of the Fermi sea drawn in 
the phase space (or more precisely of the support of the 
Wigner's function in the phase space). In equilibrium 
it is bounded by two lines = m(vo ± £po). The 

field mu(x, t) is interpreted as a wave propagating on the 
Fermi surface in the phase space. 

The time t c is roughly the time for the top of the 
wave-packet to travel a distance of its size I. Since 
the excess velocity is u, an overhang occurs at a time 
t c ~l/(5v ± %5p). At this time a curvature of the spec- 
trum becomes a dominant factor. Obviously, a smooth 
and small packet remains in the linear regime longer, but 
its life is finite. In the vicinity of the overhang the classi- 
cal Riemann equation becomes invalid. A shock must be 
regularized either by quantum corrections at A = 1, or, 
in the interacting case, by the neglected gradient terms. 

6. Dispersive regularization - the role of interaction. 
The role of dispersion changes in the presence of an inter- 
action. The latter gives higher gradient corrections to the 
Euler equation. In the linear regime they are small, but in 
the shock-wave regime, they stabilize growing gradients. 
This mechanism is called dispersive regularization, and 
is well known in the theory of non- linear waves A 
stabilization occurs when dispersive and nonlinear terms 
in J3J become of same order 

aVdp - Sp 2 . (6) 
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This condition determines the smallest scale of oscilla- 
tions occurring after a shock wave. It is Al ~ a/5p. 

This is an important result. Rising gradients are 
bounded by a scale exceeding the Fermi length k F Al^>l. 

As a consequence a semiclassical wave packet of in- 
teracting lD-fermions remains semiclassical even after 
entering a shock wave regime. Its evolution is described 
by a classical non-linear hydrodynamics. 

Condition © emphasizes the role of interaction. In 
its absence (A = 1) the gradient catastrophe is stabi- 
lized only by quantum effects and makes a semiclassical 
description non- valid at t ~ t c |23j . 

7. Chiral case A generic wave-packet will be sepa- 
rated into two parts - left and right modes moving away 
with sound velocities. The physics of shock waves is fea- 
tured in each chiral part, so we can treat them sepa- 
rately. Moreover, specifics of the Calogero model allows 
us to separate the chiral sectors exactly, by choosing ini- 
tial conditions as v = %p, so that the motion is only 
(right) chiral. Here v = v — a-^-V(logp)u is the shifted 
velocity Hermitian with respect to the inner product 
Under this condition two eqs. (|2I3|) become identical 



P- 



;-/>v(i°s(>k) = o. 

Z7T 
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8. Benjamin-Ono equation. Equation Q admits fur- 
ther simplification relevant for the physics of shock waves. 
The semiclassical nature of the initial wave packet insures 
that deviations of density from the mean value, po are 
small. Let us choose the chiral case, and a frame mov- 
ing with the velocity v s — npo- Then keeping gradient 
terms in Q of the lowest order, we linearize the disper- 
sion term pV(logp)# « V(Sp)n. As a result we obtain 
the celebrated Benjamin-Ono (BO) equation |2J] 



it + uVu + a — V 2 uh =0, u = nSp. 

4-7T 



(8) 



In hydrodynamics it describes interface waves in a deep 
stratified fluid [12|, Eqs. (121.' il) having a similar structure 
but capturing both chiral sectors, were called Benjamin- 
Ono equation on the double (DBO) in 3J. Both eqs. 
(|7I8|) are quantum. The field u in the quantum BO eq. 
is the gradient of a canonical real Bose field. 

9. Quantum Benjamin-Ono equation and FQHE. It 
is commonly accepted that edge states in FQHE are de- 
scribed by chiral bosons The quantum BO equation 
is a good candidate for an extension of the theory of 
edge states beyond the linear approximation. Indeed, 
it has all the necessary features of the physics of edge 
states: in the linear approximation it reduces to linear 
chiral bosons, the statistics of excitations is fractional, 
and it is Galilean invariant. Elsewhere we argue that 
the Galilean invariance and statistics uniquely determine 
both the nonlinear and the dispersion terms of the equa- 
tion. 



10. Semiclassical hydrodynamic equations, solitons. 
Below we study a semiclassical version of quantum hydro- 
dynamics. Due to a non- linear nature of quantum equa- 
tions, the semiclassical limit is more subtle than merely 
treating quantum operators as classical fields. In we 
argued that the semiclassical limit describes a collective 
motion of excitations (quasi-holes) with charge 1/A. This 
amounts to replacing the parameter a in (|2I3I7I8J) by the 
charge of excitation 1/A |2a] ■ To see this, one may no- 
tice that the parameter a is a charge of solitons of the 
classical non-linear field it. 

Solitons and more general multiphase solutions of BO 
are known 0, Q] . We have found general multiphase 
solutions and their Whitham modulations for more gen- 
eral DBO eas. (l2!3|l . We present these results elsewhere 
(ljj . Here we restrict ourselves to a simplified descrip- 
tion of a shock wave in the most physically relevant limit 
when BO ® holds. To this end we need to know only 
1-phase and 1-soliton solutions. 

The one-phase solution (compare to |13 |) reads 



dp- 6p=\lmd x log f 1 - 
7rA V 



v-K 
v-K* 



j-e{x,t)) 



(9) 



where 8(x, t) — k(x — Vt), k — m(v — v), V = |(u + v) 
and Sp = j:(K + ^) are the phase, the wave number, the 
velocity and the mean density. Parameters v, v and K 
are the moduli of the solution. 

The 1-soliton solution appears as a long wave limit 
k^O of © 



Sp — Sp 



1 



-Im 



1 



ttA x - Vt 



2m(V-§ Sp) 



The amplitude, 



and the width, 



(10) 



of 



2m(V~f Sp) 

the soliton are determined by the excess velocity V—^Sp. 
The charge of the soliton is 1/A. 

11. Whitham modulation. Although eqs. (|2I3I7I8(1 
are integrable, only quasi-periodic solutions enjoy ex- 
plicit formulas (l(ij . Solutions with generic initial data 
can not be found explicitly. However, in many physically 
motivated cases, the solutions are well-approximated by 
slowly modulated waves. There it is assumed that the 
moduli of a quasi-periodic solution (in the case of 1-phase 
solution the moduli are v, v, K) also depend on space 
and time but in a slow manner. They do not change 
much during a period of oscillation. If the scales of os- 
cillations and modulations are separated, the Whitham 
theory provides equations for the space-time dependence 
of the moduli [l7|. 

Shock waves is the most spectacular application of the 
Whitham theory. It has been developed in a seminal pa- 
per 0| on the example of the KdV equation. The idea 
of the method is the following. Outside the shock wave 
regime a smooth initial wave remains smooth, and one 
legitimately neglects the dispersion term, keeping only 
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non-linear terms. The equation obtained in this manner 
is the Riemann eq. JSJ. Its solution with initial condi- 
tion u(x,0) = f(x) reads u^(x,t) = f(x - u^(x,t) ■ t). 
By using the superscript we indicate that this is not the 
oscillatory part of the solution. A solution of (JSJ always 
overturns, and, typically becomes a three- valued func- 
tion in the interval X-(t) < x < x + (t) (FigQ]B), where 
the leading and trailing edges x±(t) arc determined by 

the condition d x vS ^ = oo. Let us order the branches as 

(o) . (o) . (o) 
u\ > u 2 > u 3 . 

In the three- valued region the dispersion is important. 
It replaces a non-physical overhang by modulated oscilla- 
tions. The latter are given by the Whitham's modulated 
solution. To leading order one chooses a modulated 1- 
phase solution. We call it u^\x,t) to emphasize that 
this solution has only one fast harmonic. It is glued with 
u(°) at x±(t): 4 0) (x_) = u ( °\x+) = uW(x+). 

A modulated 1-phase solution is given by the formula 
where three moduli v, v, K are smooth functions of 
space-time, and the phase 6{x, t) is replaced by a modu- 
lated phase 0(x, t) found from the Whitham equations 

VQ = m(v-v), 9 = y(u 2 -w 2 ). 

The Whitham equations for the moduli have again the 
form of Riemann equation (j^J 

v + vVv = 0, v + vVv = 0, K + KVK = 0. 

The boundary data for the moduli are chosen such that 
stops to oscillate at the gluing points x±. That hap- 
pens (i) when v — ► v at x = x+(t), there it' 1 ) — > K, and 
(ii) when v — > K, at x = X-(t), there uW — > v. The 
gluing conditions lead to an es pec ially simple result for 
the moduli of BO equation [bUllfljl : 

v = Ui(x, t), u = «2 (x,t), K = Us(x,t). 

12. Shock waves. The space-time dependence of the 
moduli approximately determines the entire evolution of 
a wave-packet. We summarize some features which are 
not sensitive to its initial shape. Let us choose a frame 
moving with the sound velocity. An initially smooth 
bump with a height Sp and a width I tends to overhang 
toward its bright side. Its dark side becomes smoother. 
At the time t c ~ l/n8p the bump starts to produce 
oscillations. The oscillations fill the growing interval 
x_(i) < x < x+(t). Further, at x > x+(i) the shape 
of the bump does not change. 

The leading edge runs with an excess velocity with 
respect to the sound. At large time x + ~ nSp ■ t. 

The trailing edge x- slows down with time. If ini- 
tial bump decays as Sp(x,0) ~ x~ n , then at large time 
X-(t) ~ i 1 /( n + 1 ). If the decay is exponential Sp(x, 0) ~ 
e~ x / r , then X-(t) ~ r\ogt. These are the simple conse- 
quences of the Riemann equation JSJ. 



The amplitude of oscillations is zero at the trailing 
edge and grows towards the leading edge. Furthermore, 
the period of the oscillations also grows. As a result, 
near the leading edge the oscillatory pattern resembles 
a collection of individual localized traveling pulses (soli- 
tons) - a soliton train. The excess velocity, V — ^Sp, 
of the leading soliton is nSp, which also determines the 
amplitude 45p. The hight of the leading pulse is four 
times higher than the initial bump. Since each pulse car- 
ries a quantized charge 1/A, the width of leading pulses, 
(AiiSSpX)^ 1 stays constant, while the distance between 
them is growing. 

15. Observation of shock waves and direct measure- 
ment of a fractional charge. Edge states of FQHE seems 
to be the best electronic system to observe quantum 
shock waves. In order to evaluate the scales one needs pa- 
rameters of the linearized theory of the edge - the sound 
velocity and the compressibility of the chiral boson. Al- 
though these quantities have never been measured di- 
rectly, one can suggest reasonable estimates [2(J ■ We es- 
timate the sound velocity on the edge as v s ~ 3 x 10 4 m/s. 
[22| An electronic density in the bulk 10 n cm~ 2 gives an 
estimate for the ID density on the edge po ~ 3 x 10 m -1 . 
It gives k — v s /po ~ 10~ 3 m 2 /s, and, at A = 3 estimates 
an effective mass m to be about 30 electron band masses 
in GaAs and the "Fermi" time tf = (^Po)" 1 to be ~ lps. 

A wave-packet of width I ~ 10 2 p " 1 ~ 1/irn and of height 
Sp^ 10 _1 po carrying about 10 electrons, develops a shock 
wave at time t c ~Ti?(/9o0f^ ~ lns - During the time t c the 
wave packet crosses a distance v s t c ~ 30pm. This scale 
is much smaller than a size of a typical heterostructure 
(~ 10 3 /im), and is still smaller than the typical ballistic 
length 50 — 100/im. Distinct solitons will start to appear 
right after the shock. Observation of the electric charge 
carried by a distinct soliton-pulse will provide a direct 
measurement of fractional charge. The full decay of this 
packet is about 10 3 times longer. 

These estimates suggest that non-linear effects (shock 
waves and fractionally charged soliton train) can be ob- 
served in the nanosecond range. Finally, we mention that 
quantum shock-waves in Bose systems have already been 
observed in trapped alkali atoms |21| . 
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Krichever, B. Spivak, A. Zabrodin, V. Goldman, and W. 
Kang. Our special thanks to I. Gruzberg, D. Gutman, 
R. Tcodorescu. PW. and E.B. were supported by NSF 
MRSEC DMR-0213745 and NSF DMR-0220198. AGA 
was supported by NSF DMR-0348358. 
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ADDENDUM. In this Addendum we provide some 
details of the quantum Double Benjamin-Ono equation 
and its classical solution. A fuller account will be pub- 
lished in 16]. 

1. Holomorphic Bose field. The seemingly compli- 
cated equations, ilL'l.'SI II) . possess integrable structure, fol- 
lowing from the integrability of the Calogero model. This 
structure becomes more transparent in terms a anti- 
Hermitian Q holomorphic Bose field p(z) and its current 
u = — i-^Vip(z) defined in the complex plane z = e l ^~ x 
Q ■ The Bose field is a Laurent series with respect to the 
unit circle ip = p+{z) + p-(z) having mode expansion 
outside of the unit circle 

( P+( Z ) = T Y] -V n Z n , tp- = V ~p n Z n + ^~\0gZ, 

A n ^— ' n 2 

n>0 n<0 

[Pn?V_ m ] — W>S nm . 

Corresponding currents are u_ = — i^V(^_, u + = 

The current is defined as a holomorphic operator hav- 
ing boundary values on the circle 

u+(x) = v-i-pH + ia—d x logp, (11) 

Z 47T 

u_(z±i0) = ^(ipHTp)- (12) 

Its jump is the proportional to the density u(x + iO) — 
u(x — iO) = —Kp(x), thus defining the modes of the neg- 
ative part of the Laurent series through the moments of 
coordinates p n — § e~ l ^ xn p{x)dx. 



Let u+{z) and u_(z) be positive and negative parts 
of the Laurent series for the current. We have u_ = 

f § p(x) z+e ]Zl ■ This field u_(z) is analytic inside 

and outside the unit disk, having a jump on its bound- 
ary equal to the density. The fact that the jump is real 
implies a reflection property 

(¥>_(*)) + = Ml/*)- (13) 

The field U-y{z) is an analytic field in some strip Ri < 
\z\ < i?2 surrounding the circle, having value u+(x) = 
v ~ i^pH + ia-^d x logp on the circle. Under the chiral 
condition it becomes u + (x) = ^(p — ip H ) + ia^S7 (log p — 
i(\og p) H )- In this case u + (z) becomes analytic inside the 
circle (i.e., R\ — in the chiral case). 

2. Quantum bilinear equation In terms of the Bose 
field the equations of quantum hydrodynamics are writ- 
ten in a compact wayQ. Let us introduce posi- 
tive/negative vertex operators r + = e~ v+ and t_ = 
e Xip - . Then eqs. H2I3J) are equivalently written in the 
form of quantum bilinear Hirota equation 

(iD t -£-D 2 x )T--T + = (14) 

with the reflection condition 

z-* 2N - l *\t-(1/z))1=t^z), (15) 

where $ is some constant. Hirota's derivatives are de- 
fined as D n J ■ g = (d x - d x ,) n } \x)g{x')\ x=x , . The bilin- 
ear equation is a bosonized version of Calogero model. Its 
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derivation from the original model and the derivation of 
the quantum hydrodynamic eqs. (121 111 from the bilinear 
equation is straightforward, but not effortless. Having a 
solution of these equations one finds the density and ve- 
locity p = ^:VlogT_, v — f^rlmV logr + . The classical 
version of l|14|). not bounded by a reflection condition, is 
known as MKP (modified Kadomtsev-Petviashvili) equa- 
tion. Here it appears in a quantum version. 

3. Reductions of MKP: Benjamin-Ono and Double 
Benjamin- Ono equations. The bilinear form is the most 
economical way to understand connections between eqs. 
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- The reflection condition l|15l) implies a reduction of 
MKP equation to the Double Benjamin-Ono equa- 
tion 

- Further requirement that t+ is an analytical opera- 
tor inside the unit disk implies the chiral conditions 
leading to eq. Q). 

- Finally further reduction = u + (l/z) + C, 
where C is a real constant, leads to the Benjamin- 
Ono equation © []]|. 

We publish a detailed study of this equation elsewhere 
[lfij . Here we present the few major facts. 

4- Semi-classical limit The classical limit is straight- 
forward in the bilinear form of the quantum equation. 
There we simply treat r + and r_ as classical fields satis- 
fying (O and (JT5J- 

Let us set u± — z f^rid x log(T±), then the classical bi- 
linear equation (|14|l reads 

u + uVu + — V 2 m = 
2m 

where u — u + + u_ and u — — u + ). 

5 Multiphase solution A general quasi-periodic (or N- 
multiphase) solution of classical MKP equation IL'l.'ill (in 
the frame where the center of mass of the system is at 
rest) is coded by 2N parameters - Vi,Vi (wave number 
ki = m(vi — Vi) and velocities Vi = + Vi) + vq of 



waves), N real phases Xi, 2N real positive amplitudes c~ 
and two zero modes K and vq. It is given by a determi- 
nant formula for the tau-functions 

( + e ^ x ^\ 

t±(x + iO, t) = e* e± det £„■ - cf — . (16) 

i,j<N y Vi-Vj J 

Here 6i(x, t) = ki(x — Xi — Vit) and 9- = —mK(x — yt), 
9 + = mvo(x — (i« + K)t), where fej = Vi — Uf, and 

Vi — — 2 V 0- 

The reflection property, l|15|) . further restricts the am- 
plitudes 

, (±), 2= ]± fvi-K^ ^f vj+K-^ . n&i(vj-Vi)(.Vj-Vi) 

where $ = ^r^i^i- The parameter $ determines mon- 
odromy of the wave function. It is the flux piercing the 
circle. Additional restrictions on Vi , t>, insure that a are 
real. The zero mode K is related to the mean density 
p = - ($ + A"). The mean velocity is v . 

One-phase solution was found in [l5| . 

The chiral case occurs when all zeros of t_ are located 
outside the unit disk of the complex plane z = e l ~ x , 
while all zeros and those of t + are inside the unit disk. 

Solutions of the BO eq. © arise at the limit Vi ^ ( K — > 0. 
This yields i|16|) with the amplitudes 01 

(c (±))2_ k J (n l -K\ ±1 l\^ l {v 3 -v l ){v ] -^) 
m 2 \u l -Kj Yltfifa-Vi^Vj-Vi)' 

6. Multisoliton solutions of the DBO equations j2Et 
appear as a long- wave limit of (|16|l . obtained by taking 
ki — > 0, while keeping Vi fixed: 

r±=e-^ e± det ((x-x l -V l t-iAf)8 lJ --^-^-) , 

where Af = ^ (^ y +K ^ v-K j ' ^ limit of this equa- 
tion (Vi — K)/K — > gives the multisoliton solution of 
the Benjamin-Ono equation = ± 2m(Vi-K) E3- 



